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O . Abstract 

Q>^ ■ /i+ -^ e+7 and yU+ -^ e+e+e^ processes are analyzed in detail with po- 

^ ' larized muons in supersymmetric grand unified theories. We first present 

I ' Dalitz plot distribution for /i"*" — * e"'"e"'"e~ decay based on effective La- 

D , grangian with general lepton-fiavor- violating couplings and define various 
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P- and T-odd asymmetries. We calculate branching ratios and asymme- 
tries in supersymmetric SU(5) and SO(IO) models taking into account 



H ' complex soft supersymmetry breaking terms. Imposing constraints from 

experimental bounds on the electron, neutron and atomic electric dipole 
moments, we find that the T-odd asymmetry for /i''" — > e^e^e~ can be 
15% in the SU(5) case. P-odd asymmetry with respect to muon polar- 
ization for /i+ — > e^7 varies from —20% to —100% for the SO(IO) model 
while it is +100% in the SU(5) case. We also show that the P-odd asym- 
metries in /i"*" -^ e~^e~^e~ and the ratio of /i"*" -^ e^e^e^ and /i^ -^ e^"^ 
branching fractions are useful to distinguish different models. 
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I Introduction 

In order to explore physics beyond the standard model (SM), rare decay experiments 
can play a complementary role to direct search for new particles at high energy 
frontier. Through forbidden or very suppressed processes within the minimal SM, we 
may be able to obtain information on interaction at the energy scale not accessible 
by collider experiments. Search for lepton flavor violation (LFV) is one of such 
windows to new physics. 

In recent years, LFV processes have received much attention because in the 
supersymmetric grand unified theory (SUSY GUT) the branching ratios for ^"^ — ^ 
e"'"7 and /i"*" -^ e^e^e" and the /i-e conversion rate in a nucleus can reach just below 
present experimental values [jll-Q. The present experimental upper bounds of these 
LFV processes are S(/i+ ^ e+7) < 1.2 x IQ-^^ §, B{^i+ -^ e+e+e") < 1.0 x IQ-^^g 
and cr{fi~Ti—>- e~Ti)/o"(|U~Ti^ capture) < 6.1 x 10^^'^ 0. It is possible that future 
experiments will improve the sensitivity by two or three orders of magnitude below 
the current bounds p, |[. 

In this paper we discuss the /i+ -^ e+7 and /i"*" -^ e+e+e~ processes in SUSY 
GUT. We focus on various asymmetries defined with the help of initial muon po- 
larization. Experimentally, polarized positive muons are available by the surface 
muon method because muons emitted from 7r"*"'s stopped at target surface are 100% 



polarized in the direction opposite to the muon momentum ||T0|- It is shown in 
Ref. |1TT| that the muon polarization is useful to suppress the background processes 
in the fi~^ -^ e"^7 search. As for the signal distribution of /i"^ -^ ^"'"7, the an- 
gular distribution with respect to the muon polarization can distinguish between 
/i"*" — > e^7 and /i+ -^ ^'rI- For /^^ ~^ e'^e^e^ , distribution in the Dalitz plot and 
various asymmetries defined with help of the muon polarization carry information 
on chirality and Lorentz structure of LFV couplings. In particular, we can define 



T-odd asymmetry which is sensitive to CP violation in LFV interactions [|T^. In 
the previous paper |T^ we pointed out that sizable T-odd asymmetry can occur in 
the SU(5) SUSY GUT when a CP violating phase is introduced in one of the soft 
SUSY breaking parameters, i.e. the universal trilinear scalar coupling constant Aq. 
The purpose of this paper is to give a model-independent framework for analyz- 
ing the fi~^ -^ e+7 and /i+ —>■ e+e+e~ processes and investigate specific features of 



the SU(5) and SO(IO) SUSY GUT focusing on the T-odd and other asymmetries. 
Detailed comparison of the T-odd asymmetry with the electron, neutron and Hg 
electric dipole moments (EDM) is also done introducing SUSY CP violating phases 
within the minimal supergravity (SUGRA) model. 

In section 2 we describe effective Lagrangian of the processes fi^ -^ e+7 and 
/i+ -^ e+e+e^. We introduce a P-odd asymmetry for fi^ -^ e+7 and two types 
of P-odd asymmetries and a T-odd asymmetry for /i+ -^ e^e^e^ . In section 3 we 
introduce the SU(5) and SO(IO) SUSY GUT and briefly review how LFV processes 
arise in these theories. In section 4 we present results of our numerical calculations. 
We calculate the branching ratios and the asymmetries in the SUSY SU(5) and 
SO (10) models taking into account complex soft SUSY breaking terms under the 
constraints imposed by the EDM experiments. We find that the T-odd asymmetry 
can be 15% in the SU(5) case while it is less than 0.01% in the SO(IO) case. We 
also show that the P-odd asymmetry for yU+ -^ e+7 varies —20% - —100% for 
the SO(IO) model and 100% for the SU(5) case. In the SU(5) case the P-odd 
asymmetries in /i+ -^ e+e+e^ can reach ±30% and the ratio of /i+ -^ e+7 and /i+ -^ 
e^e^e~ branching fractions varies over SUSY parameter space. On the contrary 
these asymmetries are smaller and the ratio of two branching fractions is almost 
constant in the SO (10) case. In Appendices, useful formulas are listed. 

II Phenomenology of the /i+ -^ e^7 and jjl^ -^ 
e+e+e~ processes 

We begin with the effective Lagrangian for /i'^ -^ e^^ and /i^ -^ e^e+e^ processes. 
Using the electromagnetic gauge invariance and the Fierz rearrangement we can 
write without loss of generality: 



^G 



+5'3(7^7''eR)(e^7^eij) + g4iWYeL)ieZ'yf,eL) 
+9^{T^l^eR){el'-ff,eL) + g%{'mi^eL){e^l^eR) + /i.e.}, (1) 



where Gp is the Fermi couphng constant and m^ is the muon mass. The chirahty 
projection is defined by the projection operators Pr = ^-^^ and Pl = ^—^- cr^^ is 
defined as a^^ = |[7^,7jy]. Al^Ai^) is the dimensionless photon-penguin couphng 
constant which contributes to /i"*" — > e^7 (/i+ -^ ^rI)- These couphngs also induce 
the /i"*" — > e~^e^e^ process, gi^s (2 = 1—6) are dimensionless four-fermion coupling 
constants which only contribute to /i"*" — > e~^e~^e~. Al^r and gi {i = 1 — 6) are 
generally complex numbers and calculated based on a particular model with LFV 
interactions. 

The differential branching ratio for ix^ — ^ e"'"7 is given by: 



dcosO 

B{fi+ -. e+7) 



1927i'{\ALr{l + cos^) + \Ar\'{1 - cos^)} (2) 

{1 + A(/i+ -^ e+7)P cos 6}, (3) 
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where the total branching ratio for /i"*" -^ e+7 {B{fi~^ -^ e"'"7)) and the P-odd 
asymmetry {A{fi~^ — > e"'"7)) are defined as 



5(^+^6+7) = 38An\\AL\'+\AR\'), (4) 



Here P is the muon polarization and 6 is the angle between the positron momentum 
and the polarization direction. 

Kinematics of the /x"*" -^ e^e^e' process with a polarized muon is determined 
by two energy variables of decay positrons and two angle variables which indicate 
the direction of the muon polarization with respect to the decay plane. In Fig. ^ we 
take the z-axis as the direction of the decay electron momentum {p^) and the z-x 
plane as the decay plane. Polar angles (6', ys) (0 < 6* < tt, < (/? < 27r) indicate the 
direction of the muon polarization P. We take a convention that the decay positron 
having larger energy is named positron 1 and the other is positron 2 and {pi)x > 0. 
We define the energy variables as xi = ^^ and xo = ^^ where Ei and Eo are the 
energy of the positron 1 and 2, respectively. In this convention [xx^x^) represents 
one point of the Dalitz plot (Fig. Q). In our calculation we neglect the electron mass 
compared to the muon mass except for the total branching ratio. In order to avoid 



logarithmic singularity we have to take into account the electron mass properly to 
evaluate the total branching ratio. 

Using the coupling constants in the Lagrangian in Eq. ([^) the differential branch- 
ing ratio for /i+ — > e^e^e~ is written as follows: 



-[Cia;i(a;i,a;2)(l + P cos9) + C2ai{xi, 0:2) (1 — Pcosi 



dxidx2d cos 6 dip 2tc 

+ C3{a2{xi,X2) + P(3i{xi,X2) cos6 + P7i(xi,X2) sin 6^ COS v?} 

+ Ci{a2{xi,X2) — PPi{xi, X2) cos 6* — P7i(a;i,a;2) sin6'cosv9} 

+ C;i{as{xi,X2) + P(32{xi,X2) COS 6' + P72(a;i,X2) sin6'cosv9} 

+ C6{Q;3(a;i,a;2) - PI32{xi,X2) cos6 - P72(a;i,a;2) sin6'cosv2} 

+ C7{Q;4(a;i,a;2)(l — P cos9) + P73(a;i, X2) sin 6' cos*/?} 

+ C8{a;4(xi,X2)(l + PcosO) — Pj3{xi,X2) sin 6* cosy)} 

+ C9{a;5(xi,X2)(l + PcosO) — Pj4^{xi,X2) sin 6* COS!/)} 

+ C'io{a5(a;i,a;2)(l - PcosO) + P'^i{xi,X2) sin 6^ cos 93} 

+ CiiP73(a;i,X2) sin6'sinv9 - Ci2P74(a;i,X2) sin6'sinv9], (6) 

where Cj are expressed by the coupling constants gi {i = 1 — Q) and A^r as: 

lo lb 

C3 = \g^\\ C4 = |^6|',C'5 = \eAR\\ Ce = leA^p, 

C7 = Re{eARgl), Cg = Pe(eAi^*), Cg = Re{eARgl), Cio = Pe(eAi^*), 
Cii = Im{eARgl + eALgl), C12 = Im{eARgl + eA^^*), (7) 

where e(> 0) is the positron charge and P is the magnitude of the polarization 
vector. Functions Oj, /?« and 7^ are defined as: 



ai{xi,X2) 

a2{Xi,X2) 
a3{xi,X2) 

ai{xi,X2) 

a5{xi,X2) 
/?l(xi,X2) 

/52(a;i,a;2) 



8(2 - Xi - X2)(Xi + X2- 1), 

2{2;i(l-Xi) +X2(1-X2)}, 

.2x^-2x2 + 1 2a;f- 2x1 + 1 



32(xi+X2-l), 



+ 



1 - X2 



}, 



(2-xi -X2) 



X 



(1 - Xi)(l - X2)(2 - Xi - X2) 

{2(xi + X2)(x^ + X2) — 4(xi + X2)(2x^ + X1X2 + 2x2) 
+ (19x^ + 30x1X2 + 19x^) - 12(2xi + 2x2 - 1)}, 



7i(a;i,a;2) 

'l2{Xl,X2) 

73(a;i,a;2) 

74iXi,X2) 



1 - Xi)(l - X2)(Xi + X2 - 1)(X2 - Xi) 



(2 -xi -X2) 



32. 



(X1+X2-I) (Xi +X2 - 1)(X2 - Xi) 
\ (1 -Xi)(l -X2) (2-X1-X2) 



16, 



(Xi + X2 - 1) 
\{l-X,){l-X2: 



■(Xi +X2 - 1)(X2 



X 



1)7 



(Xi +X2 -1) 
^(1-Xi)(l-X2) 



(2 -Xi -X2)(X2 - Xi). 



(9) 
(10) 
(11) 



K2-X1-X2), (12) 

(Xi + X2){xl + x1) - 3(xi + X2)^ + 6(xi +X2) - 4 



(13) 



(14) 
(15) 

(16) 
(17) 



In Eq. (P) there are three classes of terms: the first contribution arises from 
the four-fermion couphng constants (Ci_4) and the second from the photon-penguin 
couphng constants (C5 g) and the third from interferences between the four-fermion 
couplings and the photon-penguin couplings (C7_i2). The angular dependence with 
respect to the polarization direction is classified into four types, namely, terms pro- 



portional to (i) 1, (ii) cos^, (iii) sin ^ cosy?, and (iv) sin 6* sin 99. Under the parity 
operation (P), 9, ip transform as: 

e -^ TT-e, 

^ \ 37r-v? (7r< v?<27r) ' ^^' 

so that terms proportional to (ii) and (iii) are P-odd. On the other hand the time 
reversal operation (T) induces the following transformation: 

9^9, v? ^ 27r - V?. (20) 

Thus only terms proportional to Cu and C12 are T-odd quantities. Notice that these 
terms are given by imaginary parts of interference terms between photon-penguin 
and four-fermion coupling constants. This means that effects of CP violation can 
be seen only through a phase difference between these two coupling constants. 

It is convenient to define integrated asymmetries in order to separate four 
angular dependences, although in principle we can determine Cj separately by fitting 
experimental data in full phase space. In the Dalitz plot, ^3 and [32 have a singularity 
as ^_^ in the region near the kinematical boundary (a;i^2 ~ !)• 72, 73 and 74 
have a weaker singularity as , ^ 0:3, 82, and 72 arise as square of photon- 

penguin amplitudes whereas 73 and 74 from interferences between photon-penguin 
and four-fermion terms. On the contrary, contributions from square of the four- 
fermion coupling constants have no singularity on the edge and have a rather flat 
shape. These singular behaviors are cut off if we take into account the electron mass. 
To show this behavior explicitly, we first integrate over smaller positron energy X2 
fixing the larger positron energy xi and define the following differential branching 
ratio and three types of asymmetries ap^,ap^ and ax as a function of the larger 
positron energy Xi (| < Xi < 1): 



/ dx2 / dcos9 / dip- 

Jl-xi J-l Jo C 



dxi Ji-xi J-i Jo dxidx2d cos 9 dip 

= 3{{C^ + C2)F^{x^) + {Cs + C,)F2{x^) 

+ {C, + C,)F,{x^) + {Cr + C^)F,{x,) 



(C9 + Cio)F5(xi)}, (21) 

dB 



1 rx\ rl pItt 

- dx2 dcosO / d^——-— — -) 

Ji-xi J -I Jo dxidx2dcosud(p 



dxi 

+(C5 - Ce)G2{x,) - (Cr - Cs)F,{x,) 

+ (C9-Cio)F5(xi)}, (22) 

-1 , r^ , /■! „ rf . dB 



«^2(a;i) = pdBix,) ([ '^^dx2 [ ^d cose [^' dip- 



dx\ 

3 



dxidx2d COS 9 dip 

^1 , /"i , . /■i'^ , dB 

dxidx2dcos6dip 

^1 . /■! „ a2^ . dB 



— / dx2 / dcosO I dip 

Jl-xi J-1 -/f 

fXi rl /'27r 

+ / dx2 dcosO dip ^ , 

Ji-xi J-i J|7r dxidx2d COS Udipi 



= l^{iC,-C,)H^{x^) + iC,-Ce)H2ix,) 

dx\ 

+{C7 - Cs)Hs{x,) - (Cg - Cio)i^4(Xl)}, (23) 

^r(a;i) = — 77^7-t(/ ^^2/ dcosQ \ dip 



pdB{xi)_ J^_^^ J_^ Jq dxidx2dcosOdip 

■ix\ 

^1 . /•! „ /•2^ . dB 



pxi ri rZTV 

/ dx2 / dcosO / (iy9- 



/tt dxidx2d COS ^(ii/? ' 

3 1 



2 rfB(a'i) 



j^{Cni/3(a;i)-Ci2i^4(xi)}. (24) 



In these formulas, Fj, Gi and ifj are functions of the variable xi and their analytic 
forms are found in Appendix ^ dx ^ CLPii^i), ap^ixi) and axixi) are defined to 
extract terms (i)-(iv) with different angular dependences and axixi) is the T-odd 
quantity. In the above expression ^3(^1) in J'^^' and G2{xi) in ap-^ have j^ 



singularity. Introducing the cutoff 6 for variable xi and integrating over | < xi < 
1 — 6, we define the integrated branching ratio B and three asymmetries Ap^,Ap^ 
and At- 

"1-^ dB{xi) 



B[S] = I dxv 



dxi 
= 3{(Ci + C2)h[5] + (^3 + C,)h[5] + (Cs + C,)h[5] 

+{C7 + C,)h[5] + (Cg + C7io)/5[5]}, (25) 

1 r^-^ dB 

ApM = ^y, rfx,a,(xO^(xO 

= ^{(Ci - C2)/i[5] + (C3 - C4) Ji[5] + (C5 - C,)U5] 

-{Cj - C,)h[5] + (Cg - Cio)/5[5]}, (26) 



1 /■i-<* dB 

-{{Cs - C,)K,[S] + (C5 - C,)K2[5] + (C7 - Cs)K,[5] 



2B[5] 

-(C,-C^o)K,[6]}, (27) 

1 r^-s . , , t/5 



At[5] = ;^y. rfa;,a3(x0^(x0 

= ^{^ii^3['5]-C7i2i^4[<5]}. (28) 

Jj, Jj and Ki are functions of the cutoff 6 and their analytic forms are also found 
in Appendix ^. Note that I^IS] and J2[S] have a logarithmic singularity at 5 = 0. 
Because of this logarithmic dependence, the terms \Al\'^ and |Ai?,P dominate over 
other terms in the branching ratio if coupling constants eAi, eAp and Qi have sim- 
ilar magnitudes. On the other hand the numerator of At does not have a singular 
behavior so that At itself is suppressed when we take very small 6. In the lat- 
ter analysis of SUSY GUT cases we introduce the cutoff 6 to optimize the T-odd 
asymmetry. 

We have to take into account the electron mass properly to get precise value of 



total branching ratio. If the photon-penguin contribution dominates the branching 
ratio, we can derive a model-independent relation between the two branching ratios 



P: 










i?(/.+ - 
i?(/i+ 


-^e+7) ~ Stt^ ^ml' 
~ 0.0061, 


4^' 



(29) 

where a is the fine structure constant. Neglecting the terms suppressed by ^^, the 
total branching ratio is , therefore, given by: 

T77, 11 

i?(/i+ ^ e+e+e-) = 2(Ci + C2) + (C3 + C4) + 32{log(-|) - -}(C5 + Cg) 

+16(C7 + C8) + 8(C9 + Cio). (30) 

III SUSY GUT and LFV 

In this section we introduce SU(5) and SO(IO) SUSY GUT and discuss LFV pro- 
cesses. We assume that SUSY is broken explicitly at the Planck scale with soft 
SUSY breaking terms and that these terms have universal structure with respect to 
the fiavor indices as suggested by the minimal SUGRA model. First, we discuss the 
LFV process in the SU(5) SUSY GUT and introduce the SO(IO) SUSY GUT in the 
next subsection. 

III.l SU(5) SUSY GUT 

In the SU(5) SUSY GUT, we have three generations of 10(T) and 5(F) representa- 
tions of SU(5) as matter fields and ^{H) and 5(if) representations of Higgs fields. 
The Yukawa superpotential and the soft SUSY breaking Lagrangian are written as 
follows: 



>V5^(5) = \{yu)^,T{^^H + {yd),^F,T,H, (31) 






-{^{Au)ijT,TjH + mo{Ad)ijFiT,H + -M^X^^^l + h.c.}. (32) 

o / 

where i, j are generation indices. T, F are scalar components of the superfields T, 
F. 

At the Planck scale these soft SUSY breaking parameters satisfy flavor-blind 
universal conditions which are implied in the minimal SUGRA model: 



2 2 2i 2 2 2 

rrij, = Trip = mQl, rrijj = m-^ = mg, 

{Au)ij = Ao{yu)ij, {Ad)ij = AQ{yd)ij. (33) 

With these conditions the lepton and slepton mass matrices can be diagonalized 
simultaneously at the Planck scale, and therefore there is no LFV at this scale. 
However, these conditions receive corrections from the renormalization effect be- 
tween the Planck scale and the GUT scale mainly due to the large top Yukawa 
coupling constant. As a result the magnitude of the 3-3 element of the mass matrix 
for 10 scalar fields becomes smaller than 1-1 and 2-2 elements. In the basis where 
yu is diagonalized at the Plank scale, the mass matrix for the 10 scalar fields at the 
GUT scale is approximately given by: 



m^ ~ I m^ 






57r^ Mr 



Am^ c^ -^l(z/«)33lX(3 + |^r)ln(^), (34) 



where Mp and Mq denote the reduced Planck mass (~ 2 x 10^^GeV)and the GUT 
scale (~ 2 X lO^^GeV). This correction amounts to about 50% of their original values 
and the lepton and slepton mass matrices are no longer diagonalized simultaneously. 
This becomes a source of LFV which could induce observable effects in /i+ -^ e+7[|ll. 
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The SU(5) symmetry is broken to the SM groups at the GUT scale, and after 
integrating out heavy fields the effective theory becomes the minimal supersym- 
metric standard model (MSSM). The superpotential and the soft SUSY breaking 
Lagrangian for the MSSM are written as follows: 

+e''^ym,H2jJtQjp + e^^^iH.^H,^, (35) 

^soft = -{ml)ijE*Ej - {ml)ijL*Lj - {ml)ijD*Dj 



+-M,BnBL + -M^WrWl + -M^GrGl + h.c.]. (36) 

In this formula e"^ is defined as e^^ = e^^ = 0, e^^ = — e^^ = 1. At the GUT scale 
these parameters satisfy the GUT relations: 

Ve = Vd, (37) 

Ae = Al (38) 

2T IT 1111 1 

m^ = Tjiu = tUq = rrij,, m^ = rrij^ = m-^, 

mjj^ = m|, mj,^ = m\, M^ = M., = M^ = M5. (39) 

In the basis where y^ is diagonalized at the Planck scale, yu at the GUT scale still 
approximately remains diagonal. In this basis, y^ is diagonalized in the following 

way: 



yny^yl = diagonal, (40) 

where Vl and Vr are unitary matrices and using Eq. (|37|) Vr is given by: 

(VRh = (V^mU (41) 
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where V^^v/ is the Kobayashi-Maskawa (KM) matrix at the GUT scale. 

It is useful to make unitary transformations on Ei and Lj to go to the basis 
where i/e is diagonalized at the GUT scale. In the new basis the off-diagonal element 
of m^ is given by: 

«h ^ -^,{VhMVLf)y{yuWml{3 + \Ao\') In(^). (42) 

The off-diagonal element of the slepton mass matrix becomes a source of LFV. 

In the actual numerical analysis, we solved the MSSM renormalization group 
equation from the GUT scale to the electroweak scale and determine the masses and 
mixings for SUSY particles. We also require the electroweak symmetry breaking 
occurs properly to give the correct Z-boson mass. From the MSSM Lagrangian 
at the electroweak scale we can derive the LFV coupling constants A^^Ft and Qi^q 
through 1-loop diagrams involving slepton, gaugino and higgsino. The complete 
formulas are given in Appendix p.2| . 

In the SU(5) model, only the right-handed slepton mass matrix can develop off- 
diagonal terms if the ratio of vacuum expectation values of two Higgs fields (tan j3 = 
jj^) is not very large. In such a case only A^, g^ and g^ have sizable contributions. 
Restricting to small or moderate tan/5 cases, all effective coupling constants are 
proportional to the product of the KM matrix element A,- = (^^M)32(^xAf )3i since 
the LFV transition occurs through (m|;)2i or (^|;)32("^l;)3i- This situation does not 
change even if we take into account the LFV transition due to the left-right mixing 
of the slepton mass matrix. This means that the CP violating phase of Yukawa 
coupling constants cannot make a phase difference between Al and g^, or Al and 
g^, and therefore the T-odd asymmetry A^ cannot appear from this source. 

There is another important source of CP violating phases in soft SUSY break- 
ing terms. Within the SUGRA model, we can introduce four phases: phases of M5, 
Aq, B and yU, but not all of them are physically independent. By field redefinition, 
we can take the phases of Aq and /i as independent phases. If we take into account 
these phases. At can be generated. Since these phases also induce the electron, neu- 
tron and Hg EDMsfl^, 0, we take into account these EDM constraints to obtain 
allowed region of SUSY phases. 

Up to now we consider that the Yukawa coupling constants are given by Eq. 
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(pip, so that the lepton and down- type quark Yukawa couphng constants are related 
at the GUT scale by Eq. (|37D. On the other hand, it is known that this relation does 
not reproduce realistic mass relations for charged leptons and down-type quarks 
in the first and second generations. It is therefore important to study how the 
prediction for LFV processes depends on details of the origin of the Yukawa coupling 
constant in the MSSM Lagrangian. One way to generate realistic mass matrix is 
to introduce higher dimensional operators in the SU(5) superpotential. Once this 
is done the simple relationship between the charged lepton and down-type quark 
Yukawa coupling constants does not hold. Although the effect of higher dimensional 
operators is suppressed by O (|f^), masses and mixings for the first and second 
generations can receive large corrections to the GUT relation. If tan/5 is not very 



large, LFV is still induced only for the right-handed slepton sector and Eq. ( p2D 
holds with replacement of V^^f tiy V^ which is not necessary related to the KM 
matrix elements. In the followings, therefore, we treat At as a free parameter. Since 
the fi~^ — >■ e+7 and the /i"*" — * e^e'^e' branching ratios are proportional to jA^P, we 
present these branching ratios divided by |At-P. If tan /5 is as large as 30, the bottom 
Yukawa coupling constant can induce the LFV in the left-handed slepton sector. In 
such a case, if we include the effect of higher dimensional operators at the GUT 
scale, there are photon-penguin diagrams which are proportional to 771^ and these 
contributions tend to dominate over other contributions as shown in [l^. Because 



the LFV branching ratios depend on many unknown parameters in such a case, we 
do not consider this possibility here. 

III.2 SO(IO) SUSY GUT 

In the minimal SO(IO) model, we assume three generations of 16 representation 
matter fields (\I'i) and two 10 representation Higgs fields ($u,$d) of SO(IO). The 
Yukawa superpotential and the soft SUSY breaking Lagrangian are written as fol- 
lows: 



Wso(io) = -{yuh^i^u^j + ^iVdh'^i^d^j, (43) 



Csoft = -(m|),,^l^,-m|^<l>t$„-m|^$|i<I> 



''■■d 
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^ ^0 / . N T ^ T "^0 / ^ X T ^ T 1 



+h.c.}. (44) 

At the Planck scale, we have the universal boundary conditions: 

ml = mil, m%^ = m|^ = ml, {Au)ij = Ao{yu)ij, {Ad)ij = Ao{yd)ij. (45) 

In contrast with the SU(5) SUSY GUT, all matter fields are unified in a single 
representation \1' of SO(IO) and masses of all squarks and sleptons of the third 
generation receive a large correction due to the renormalization effect by the top 
Yukawa coupling constant. In the ?/u-diagonalized basis, difference between the 
mass of the third generation sfermion and that of the first and second generation is 
given by: 

Ami ^ -^,\{yu)33\'ml{3 + \Ao\')\nA. (46) 

At the GUT scale, the initial conditions for the parameters of MSSM Lagrangian in 
Eqs. ( P^D and ( P^ ) at the GUT scale are written as follows: 

Ve = Vd, (47) 

Ae = Ad, (48) 

ml = ml = ml) = mlj = ml = m|, 

m|^ = m|^, ml^ = m|^, M, = M^ = M^ = M^o, (49) 

where the symmetric matrix ye can be expressed as: 

ye = U^PyeU, 



P = e'^^ , (50) 
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where jje is a real diagonal matrix, and therefore the unitary matrix U is related to 
the KM matrix at the GUT scale as: 

U = V^^^. (51) 

If we go to the j/e-diagonalized basis at the GUT scale, the off-diagonal elements of 
slepton mass matrices become as follows: 



«)., ^ -^e-^('^'-^^n^^M)3.(V^^M)3,l(2/033pm^(3+|Aonin(^), (52) 

OTT J-Vic 



'G 



A\, ^ -^(^^M);.(^^M)3,l(l/.)33pm^(3 + |AoP)ln(^; 



Since the left-handed slepton also has the LFV effect in the case of the SO (10) SUSY 
GUT, there are dominant photon-penguin diagrams which are proportional to m,- 
in the slepton left-right mixing as discussed in p[ . 



In addition to the KM phase, there are two physical phases in Eq. (50) up to a 



overall phase. A combination of these phases and the KM phase is responsible to the 



electron EDM[18|,g]. If the photon-penguin diagram proportional to 771^ dominates 
in the /i^ -^ e+7 amplitude, there is a simple relation between the electron EDM 
and the /i"*" -^ e"'"7 branching ratio 0. Defining a phase as 

/m[e^(^3-0.){(y^^^)3^(y^^^)*j2] ^ |(y^^,)3,(Vl^)*3|2sin0, (54) 

the relation is given by: 



141 = 1-3V ,^_,o I sin 01 (lO-^^e- cm). (55) 



Later we see that the diagram proportional to m^ does not necessarily dominate 
over other diagrams. In such a case the above relation does not hold. 

IV Results of numerical calculations 

We present results of our numerical analysis on //+ -^ 6^7 and /i+ -^ e+e+e^ 
processes for the SU(5) and SO(IO) SUSY GUT. We also calculate the electron, 
neutron and Hg EDMs as constraints on the CP violating phases of the soft SUSY 
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breaking terms. Following the procedure discussed in the previous section, we solve 
the renormalization group equations with the universal condition for the SUSY 
breaking terms at the Plank scale. Though the approximate formulas for the slepton 
mass difference are given in the previous section to explain qualitative features, we 
solve the renormalization group equations from the Planck scale to the electroweak 
scale numerically taking into account the full flavor-mixing matrix for fermions and 
sfermions. To determine allowed range of SUSY parameter space we use the results 
of various SUSY particle searches at LEP and Tevatron and the branching ratio 
B{b —>■ s'-f). The details on these constraints are described in Ref. [|1910- We take 



the top quark mass as rrit = 175 GeV. Because we calculate the LFV branching ratios 
divided by |A^p, the result is almost independent of the KM matrix elements. For 
definiteness, we use the input parameters of the KM matrix elements as |(V/fAf)cfe| = 
0.041, |(Vft-M)td| = 0.006 and |(Vft-M)us| = 0.22. Requiring the radiative electroweak 
symmetry breaking the free parameters of the SUGRA model can be taken as tan/3, 
M2, mo, \Aq\ and the phase of Aq {9ao) and that of /x (6*^). 

IV.l SU(5) GUT 

Let us first discuss the case without the CP violating phases in the SU(5) GUT. In 
Fig. ^ we present the following quantities 

B{fi+ -^ e+7) B{fi+ -^ e+e+e-) B{fi+ -^ e+e+e-) 
\Xj^ ' \XJ^ ' B{fi+^e+-f) ' 

A(/x+^e+7), Ap, ,Ap„ (56) 

in the plane of mg^ and |y4o| for tan/3 = 3, M2 = 150GeV, 9ao = 9^ = 0. Here A^ 
is defined by the mixing matrix which diagonalizes the right-handed slepton mass 
matrix at the electroweak scale in the basis where the charged lepton mass matrix 
is diagonalized. For the asymmetries we take the cutoff parameter 6 = 0.02. If 
|A^| = 10"^, B{fi^ -^ e+7) can be 10~^^ and B{fi -^ e+e+e~) can be 10~^^ level, but 
if Xr is given by the corresponding KM matrix element, \Xr\ becomes (3 - 5) x 10""^, 
so that the branching ratios are smaller by three orders of magnitude. In Fig. ^(c) 
the ratio of two branching fractions is shown. If the photon-penguin contribution 
dominates over four-fermion ones this ratio is given by Eq. (pOj). We can see that 



§The branching ratio B{b -^ 57) is updated as 2.0 x 10"-* < B{B -^ X^-f) < 4.5 x 10"'' [M 
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for large parameter region the ratio is enhanced. In particular, near mg^ = 400-600 
GeV almost exact cancellation occurs for the photon-penguin amplitudes 0. In 
Fig. |^(d) A{fi^ -^ e+7) is shown. It is close to 100% except for small region where 
the almost exact cancellation occurs. The P-odd asymmetries Ap-^ and Ap^ are 
shown in Fig. |^(e) and Fig. 0(f). Ap-^ changes from —30% to 40% and Ap^ changes 
from —10% to 15%. For 6 = 0.02 the asymmetries Ap-^ and Ap^ are expressed as: 



Ap, ~ ;^{0.6(Ci-C2)-0.12(C3-C4)+5.6(C5-C6) 

ZJd 

-4.7(^7-^8) + 2.5(^9 -Cio)}, (57) 

Ap, ~ A{o.i(C3-C4) + 10(C5-C6) 

+2{C7-Cs)-1.6{C,-C,o)}. (58) 

In the SU(5) case, because only g^, g^ and Al have sizable contributions, we obtain 
the following expressions: 



^{0.6|(?3|'-0.16|55r-5.6|eAi|' 

+A.7Re{eALgl) - 2.5i?e(e/li^*)}, (59) 

Ap, ^ ^{0.1|(75|'-10|eAi|2 

-2Re{eALgl) + l.QRe{eALgl)}. (60) 

In the above formulas we can see that the coefficients for I^lP, Re{ALgD and 
Rei^Aigl) are large. Therefore these asymmetries represent the dependence of square 
of photon-penguin terms and interference terms. It is interesting to see that we 
can over-determine the three coupling constants (73, g^ and Al from observables 
B{^+ -^ e+7), B{^+ -^ e+e+e"), Ap^ and Ap^ if we assume the SU(5) SUSY GUT 
without the SUSY CP violating phases. For example, we can determine (73, g^ and 
Al from the three observables -B(/i^ — > ^^l)i B{^ -^ e~^e^e~) and Ap^, then, Ap^ 
can be predicted. In addition we should have A(/i+ — > e+7) = 100% and At = 0. 
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Next, we include the SUSY CP violating phases and discuss EDM constraints 
and T-odd asymmetry. We calculate the electron and neutron EDMs according to 
Ref. ||2T[. Discussion on QCD correction is given in Appendix y. For the Hg EDM, 



we use the result of Ref. [|l6l . dng is given: 

dHg = -iC^' - Cf - 0.012Cf ) X 3.2 ■ lO-^e, (61) 

where C^ , C^ and Cf are chromomagnetic moments discussed in Appendix 0. 

In order to see 9ao and 6'^ dependences on the EDMs and At, we first show 
these quantities for a specific set of SUSY parameters. In Fig. ||, the electron, 
neutron and Hg EDMs and At are shown for tan/? = 3, M2 = 300GeV, TTig^ = 
650GeV, l^ol = 1 in the parameter region —n < Oaq < tt and — O.OSvr < 6^ < O.OSvr. 
The experimental bounds on the EDMs are given by \df.\ < 4 x 10^^^(e ■ cm) p2| , 
\dn\ < 0.63 X 10-25(e • cm) ^ and \dHg\ < 9 x 10-28(e ■ cm) |2§. As is well known 
in Ref. |^ the EDMs are very sensitive to 6^, so that 6^ is strongly constrained. 



On the other hand 6ao can be large. In this particular parameter set, 6^0 = f is 
not excluded by three EDM constraints. Maximum value of the T-odd asymmetry 
At in allowed region in this figure is 15%. Note that At is proportional to sin^^g 
in a good approximation because the magnitude of 6^ is strongly constraint by the 
EDMs. 

In Fig. 1^ we show the quantities in Eq.(|56|) and At for tan/5 = 3, M2 = 
300GeV, 9ao = |, ^^ = 0. We also show the constraints from the electron, neutron 
and Hg EDMs. Within the EDM constraints At can be 10%. As discussed in Fig. ^, 
when we vary ^^ around 6^ = 0, the EDM values change considerably but At is 
almost constant. Therefore the allowed region by the EDM constraints moves in 
the Fig. |] if we take 6^ as slightly different value from 0. On the other hand the 
contours for branching ratios and the asymmetries in this figure are almost exactly 
the same. In this figure we also show the parameter region which is not allowed 
by the EDM constraints even if we change 6^ around 6'^ = for 6ao = f • Within 
the allowed region, the maximum value of At is 15%. A similar plots are shown for 
tan/5 = 10 in Fig. ^ In this case also the maximum value of At is about 15%. Note 
that, in the case with the CP violating phases, we can still determine the complex 
coupling constants g^, g^ and A^ up to a total phase from the two branching ratios 
B{fi~^ — > e~^-y), B{fi~^ -^ e+e"'"e~) and three asymmetries Ap^, Ap^ and At- 
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IV.2 SO(IO) GUT 

In the SO(IO) case, from Eq. (|50|) , there are two physical phases which contribute 
to the EDMs and /i"*" -^ e"^7 amphtudes. In the /i"*" -^ e"^7 amplitudes the term pro- 
portional to rrir has a dependence of e*(''^3-</'2)(-y^^)g2{(V^jv,f)33}^(^^Af)3i and other 
contributions are proportional to (V^A/)32(^xAf)3i- Therefore, the branching ratio 
/i"*" — >■ e+7 depends on the relative phase of two terms. In the foUowings we consider 
the case where there is no relative phase so that the amplitude is proportional to 
A^. Also we do not consider EDM constraints from Eq. (|55|) explicitly since this 
can be suppressed when (p is small. 

In Fig. the branching ratios and the asymmetries are shown for the SO (10) 
model. We first show the case without the SUSY CP violating phases. Input SUSY 
parameters are taken as tan/5 = 3, M2 = ISOGeV, 9ao = and ^^ = 0. We see that 
B{fi~^ -^ e~^'y)/\Xr\'^ can be 10~^. This value is enhanced by 2-4 orders of magnitude 
compared to the SU(5) case. The ratio of two branching fractions is almost constant 
because the photon-penguin diagrams give dominant contributions to /i"*" — > e~^e~^e~ . 
The /i"*" — i> e"'"7 asymmetry A{fi~^ -^ ^~^l) varies from —20% to —90%. This is in 
contrast to the previous belief that Al and Aji have a similar magnitude in this 
model. Although the diagram proportional to rrir gives the same contribution to 
the Al and Aji, there is a chargino loop diagram which only contributes to A^i. 
In spite of no rrir enhancement, the contribution from the latter diagram can be 
comparable to that from the former one, especially when the slepton mass is larger 
than the chargino mass. In Fig. |^(e), (f) the P-odd asymmetries for /i"*" -^ e~^e~^e~ 
are shown and these asymmetries are small compared to the SU(5) case. Api is less 
than 10% and Ap2 is less than 14%. In this case C5 and Cq terms dominate in Eqs. 
(0) and (^) so that these asymmetries are proportional to A{fi~^ —^ ^"'"7). We 
have also investigated the case with tan/5 = 10. We found the parity asymmetry for 
yU+ — > e+7 and Api, Ap2 have a similar magnitudes as Fig. 0, namely, v4(/i+ — > e+7) 

varies -20% 100%, Api varies 2% - 10% and Ap2 varies 4% - 16% in the same 

parameter space. 

In Fig. Ijwe consider the case with the SUSY CP violating phase and take input 
parameters as tan/5 = 3, M2 = 300GeV, 9ao = f and 9^ = 0. The branching ratio 
and other asymmetries have similar magnitudes compared to the case in Fig. |^. We 
can see that the T-odd asymmetry At is less than 0.01% because only the photon- 



19 



penguin amplitude becomes large. 
Some remarks are in order: 

1. When we take into account the phase in Eq. (|5^, the EDM is generated as 
discussed in Eq. (|55|) . We note that the T-odd asymmetry cannot be large 
even in such a case because the photon-penguin diagram dominates over the 
four-fermion contributions. 

2. If the /i^ -^ e^'y asymmetry is sizable, the simple relationship between the 
EDM and the /i"^ —>■ e"'"7 branching ratio as in Eq. (|55|) does not hold. This 
is because the EDM amplitude is no longer proportional to the /i"*" -^ e"'"7 
amplitude due to the chargino loop contribution. 

3. Even if we include relative phases between the term proportional to ttIt- and 
other contribution in the /i"*" — > e"'"7 amplitude, we expect large A{fi'^ —>■ e+7) 
as long as two contributions have a similar magnitude. We have numerically 
checked that the asymmetry varies 0% 100% by changing the relative phase. 

IV. 3 Differential branching ratio and asymmetries 

Up to now we only discussed the integrated branching ratio and asymmetries of 
fi~^ —* e^e^e^ . In the actual experiment, the differential quantities are useful to 



distinguish different models. For example, in Fig. ^ and |T0|, we show the differential 
branching ratio and asymmetries for a particular parameter set in the SU(5) and 
SO(IO) models. |^, ap^, ap^ and a^ are plotted for the parameter set of tan/5 = 3, 
mg^ = 700 GeV, M2 = 300 GeV, \Ao\ = 0.5, Oao = f and 6^, = 0. We can see clear 
difference between the SU(5) and SO(IO) models. The differential branching has the 
steep peak near xi = 1 for the SO(IO) case whereas the distribution is broader for 
the SU(5) case. This is because the photon-penguin contribution has j^ behavior 
near xi = 1 and the four-fermion operators give a broad spectrum. We also see the 
T-odd asymmetry has the peak at xi close to 1. This fact arises from the ,-^^_^ 
behavior in the 73 and 74 near xi = 1. Because of this feature of distribution, we 
have chosen S = 0.02 to optimize the T-odd asymmetry. 
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V Conclusion 

We developed the model-independent formalism for the process yU"*" —>■ e+7 and 
yU"*" — > e~^e~^e~ with polarized muon and defined convenient observables such as the 
P-odd and T-odd asymmetries. Using explicit calculation based on the SU(5) and 
SO(IO) SUSY GUT, we show that various combination of LFV coupling constants 
can be determined from the measurement of branching ratio and asymmetries. In the 
SO(IO) case the P-odd asymmetry in fi'^ — > e"^7 varies from —20% to —100% whereas 
it is +100% for the SU(5) case. With the branching ratio and the P-odd asymmetries 
in the fi^ -^ e^e^e~ process, we can over-determine the coupling constants in the 
effective Lagrangian in the SU(5) SUSY GUT if there is no SUSY CP violating 
phases. We also calculated the T-odd asymmetry in the yU"*" — > e^e+e^ process with 
the SUSY CP violating phases and compare it with the neutron, electron and Hg 
EDMs. In the SU(5) case we can still determine these coupling constants using 
additional information of the T-odd asymmetry. The T-odd asymmetry can reach 
15% within the constraints of the EDMs. In the SO(IO) case the T-odd asymmetry 
is small as a result of the dominance of photon-penguin diagram. These results are 
summarized in Table ||. We stress that although the magnitude of the branching 
ratio has a large uncertainty due to the unknown parameter A^, asymmetries and 
the ratio of two branching ratios are independent of this ambiguity. Thus these 
quantities are useful to distinguish different models. 

The experimental prospects for measuring these quantities depend on the 
branching ratio. For the SO (10) model we expect the /i"*" — >■ e^7 branching ra- 
tio can be 10~^^ when the A^ is given by the corresponding KM matrix elements. 
In such a case the /i"*" -^ 6^7 asymmetry can be measurable in an experiment with 
a sensitivity of 10~^^ level. For the SU(5) model, to get the /i"*" -^ e^^ branching 
ratio of order 10^^^ and yU+ — > e^e^e~ branching ratio of 10"^'', we have to assume 
At- is larger than a several times 10"'^. If the branching ratio turns out to be as 
large, the ^i^ -^ e+e+e^ experiments with a sensitivity of 10~^^ level could reveal 
various asymmetries. Because various asymmetries are defined with respect to muon 
polarization, experimental searches for [i^ -^ e"^7 and /i+ -^ e+e+e~ with polarized 
muons are very important to uncover the nature of the LFV interactions. 

The authors would like to thank Y. Kuno and J. Hisano for useful discussion. 
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SU(5) SUSY GUT 


SO(IO) SUSY GUT 


A(^+ ^ e+7) 


+100% 


-20% - -100% 


B{^l+^e+e+e-) 
B(iM+^e+-y) 


0.007-0(1) 


constant (~ 0.0062) 


Ap, 


-30% - +40% 


<10% 


Ap, 


-20% - +20% 


<15% 


At\ 


< 15% 


< 0.01% 



Table 1: Summary of the results 
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A Branching Ratio and Asymmetries 

In this appendix, we give kinematic functions which appear in the calculation of 
branching ratio and asymmetries. 





Fi{x) = 2 / dx2ai{x,X2), 


(62) 




Gi[x) = 2 / dx2^i[x,X2), 

J\-x 


(63) 




Hi{x) = ~2 dx2-fi{x,X2), 

Jl-x 


(64) 


F,{x) = 


-- -^(4a;-5)(2x-l)^ 


(65) 


F2{x) = 


= -^(2a;-l)(8a;2-8a;-l), 


(66) 


Fsix) = 


= 16 In^ )(2x^-2x + l) + -^- '-^ '-, 

1 — X 3 1 — X 


(67) 


F,{x) = 


= 32(2a;-l)^ 


(68) 


F,ix) = 


= -8(2a;-l)(2a;-3), 


(69) 


Gi(x) = 


= -16(l-a;)2ln2(l-a;) - - (2a; - l)(8a;2 - 32a; + 23), 


(70) 


G2{x) = 


= -16(2a;2 - 2a; - 7) In 2(1 - a;) + 16(2x2 - 2x + 1) In 2x 
^ 32 (2a;- l)(x2 - 13x + 13) 


(7^) 



1 — a; 



Hi{x) = 2(6-5a;)(l-a;)V2a;- 1 

.2 -3a;. 



-(7a;2 — 24a; + 16) Vl — a;arccos( 



X 



+16(l-a;)2arccos( -)}, (72) 

X 



H2{x) = -16(6-a;)V2x-l 
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5a;2 + 8a;-16 2 - 3a; /^-^a f7n\ 

—8 ^ arccosf ) — 128arccos( ), (73) 

V 1 — a; X X 

Hsix) = -\^j2^^\(XIx^ - 24a; + 4) + 2 ^^^~^^^ arccos(^^^), (74) 
3 \J\ — X X 

2 
//^(a;) = +-V2a;-l(17a;2 - 30a; + 16) 
o 

(7a;2 - 16a; + 8)a; ,2 -3a;, 

■ arccos( ). (75) 



Vl — a; a; 



/,[5] = r 'rfxF,(a;), (76) 

/•l-<5 

<^j[^] = / dxGi{x)dx, (77) 

i^if^] = / dxHi{x)dx, (78) 

•^5 



/i[<5] = ^(l + 2<5)(l-25)^ (79) 

/2[<5] = ^(1 + 2(5 -2(52)(1- 25)2, (80) 

h[6] = ^{l-6)i2-5 + 25')H^)-^il-26)il3-A5 + A5'), (81) 

16 

h[6] = y(l-25)^ (82) 

h[5] = ^(l + 5)(l-2<5)^ (83) 

Ji[5] = -1-^6 + 66' + f{\n26 -^)6' -l6\ (84) 

16 16 

J2[5] = {2 + 216 + 36'-25^)\n25+—il-5)i2-6 + 25')ln2{l-6) 

o o 

--(l-25)(49 + 685 + 4(52), (85) 
9 
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Ki[5\ = -—(8 + 86- 936^ - 2256^)Vl - 26 

o 10 

-^(5t(l -Q6- 36^) arccos(^^-^) 
3 1 — 

-—6^ B.Tccos{- -), (86) 

O 1 — 

QQ or T 

K2\6] = — {A + 96 + 6^)Vi^^-IQV6(3 + Q6- 6'') avccosi—^) 
5 1 — 

r 

+1285arccos(- -), (87) 

1 — 



Ks[6] = ^^1325(48 - 575 - 685^ + 855^) 

-LUO 

-4(1 - 6fV6 aiccosi ~, ), 

1 — 



K4[6] = -—VT^^{64 - 416 + 266^ - 856^) 

L\JO 

-2(1-6-6^ + 6^) V6aTccos{^^^). (89) 

1-0 

B LFV effective coupling constants in MSSM 

B.l MSSM Lagrangian 

We first fix our notations of tlie MSSM for tlie numerical calculation. Using 

(90) 



(91) 



(92) 
The neutralino and cliargino mass matrices are written as follows: 







V ~- 


^^2{{mf + {H 


'l)\ 


and 






-^HS' 




the charged lepton 


mass 


matrix 


is given: 








(m, 


'Jtj [Vejij p- 


cos/3. 
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c 



N 



' Br W,r H^ 



R 



^2R 



M 



N 



( B, \ 

WZL 

\ HL J 



+ h.c, 



M 



N 



I 



\ 





M2 



Ml 


—m^ sin 6*1^ cos /5 m^ cos ^vi^ cos /3 
rriz sin 6'vi/ sin /3 — n ■ , 



—niz sin 9w cos P rUz sin 6w sin /9 \ 

niz cos 6'i4/ cos /? —THz cos 6'vK sin /3 

-/i 

m^ cos 6*11/ sin {3 — /i / 



£ 



c 



( Py^ 



R H:,^)Mcy~'l^^+h.c., 



M 



c 



M2 \/2mw cos P 

^/2mw sin P n 



They are diagonalized with unitary matrices as follows: 



(93) 



(94) 



OnMnOJi = diag{m~o, m~o, m~o,m~o) 



Xi' "X2' X3' Xi^ 



Onr 

Onl 

it 



On, 



On*, 



OcrMcO'^cl = diag{m--,m^~ 
The slepton mass matrices are written as follows: 



(95) 

(96) 
(97) 
(98) 



C, 



m= 



L-t Et 



m- 



L- 

E 



m\ + m\me + ml cos2/9(— | + sin^6'vK) -;^ cos/3(moAe + ye/U* tan /?)"'' 
-^ cos P{moAe + i/efj.* tan /3) m|; + mem\ — ml cos 2/3 sin^ 6w 



(99) 
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ml = ml- -ml cos 2/3. (100) 

They are diagonalized with unitary matrices as follows: 

Ue-mjUl = diag{ml^,ml^,ml^,ml^,ml^,mlj, (101) 

U^mlUl = diag{m\,m\,m%). (102) 

The neutralino and chargino vertices for leptons and sleptons are written as follows: 

^ = el{Nl'AxPL + N^AxPR)x'Aex 
+el{Ct,xPL + C^AxPR)XAi>x 
+h.c., (103) 



Ntix = -9{V2t^n9wiO^L)Mmh+3+ J"^''"' (O^L):43(^e)x,}(104) 

VZrriw COS p 



\/2mw cos P 



CtAx = 9 /K ' ^ {OcLTA^iU^Txr (106) 

V 2m vi/ cos p 

C^Ax = ~9iOcRrAiiU.)h, (107) 
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B.2 LFV effective coupling constants 

The formulas of effective coupling constants for fi ^ e'j and yU. — »• 3e processes 
written in the MSSM variables are given in Ref.[^. We present these formulas for 
completeness with taking care of the CP violating phases. 

Each coupling constant is divided into a neutralino-charged-slepton-loop con- 
tribution and chargino-sneutrino-loop contribution. The four-fermi coupling con- 
stants are given as follows: 

g^ = g? + g':{^ = ^-Q)■ (108) 

The coupling constant gi comes only from box diagrams. 



m^o m^o do{m^o '^,m^o '^ , rriex'^ , rrigY'^) , (109) 

/9 2 3 

yi — ~ dA in Z^ Z^ ^2AX^\AY^\BY^\BX 

'^Xa'^Xb^^^'^Xa'^ ' "^Xb' ^^^^-^ ^^y'^)- (110) 

The coupling constant (73 is divided into three parts. (731 is a contribution of box 
diagrams and (732 is that of Z-penguin diagrams. (733 is a contribution of a off-shell 
photon-penguin diagrams. 



6-3 = fi-Sl + 5-32 + 5^33, (111) 






^?l = - (,A^2n_ E E {^2AX^xlY^XBY^lBxdMx\^^X%^^ex\^eY'') 



A,B=ixy=\ 



+ 2^2AxNiAYNiBYNiBx1T^x'i^X%Mm^<i^, Vfl^^^ , m^J, ^eV^)}, (112) 

^3^2 = -T7^^J^[E Y.N2AxN^Bx{^{Y^-^ABC2{m,,\m^.^\m£) 

^^^ A,B=1X=1 
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-2m^o m^o {Y^o)abCo {m-e^ ^ , m^o ^ m^o 2) } 



+ E E A^2ViVl':4V(^eJxyC2Kc.^me/,mev^)], (113) 

A=i xy=\ 

^3-3 = -TTT^-IT^ E E ^2'4xiV/:lx-^&o(-^), (114) 



11527r2G'F^iir'i 



ex ex 



/^ 2 3 

"31 ~ a A in 2-^ 2-^ ^2AX^\AY^\BY^\BX 

b47r LxF A,B=\xy=\ 

C?2('T^^-^,m^-^m^x^"^^>y^)^ (115) 



1 2 3 

^32 = -T7—2^R E EC2VCl'Sx{4(V'^-)Ai?C2(m,^^m 2,m^^ 

-■-"^ A,B=1X=1 



2n 



-2mJ-m--(y--)ABCo(mp^^rn--^m--2)}, (116) 



''Xa"''Xs^"Xh^^^"UV"'';^x j'-x^ ,"''X^ 



,2 



^/Op2 2 3 1 mt_ 

^33 = -T^Bl?rllllCkAxC{l^^h^-^). (117) 



11527r2G'i.^^^^ 



rrtr, m„ 



The coupling constant (75 is divided into three parts. (751 is a contribution of box 
diagrams and (752 is that of Z-penguin diagrams. (753 is a contribution of a off-shell 
photon-penguin diagrams. 

95 = 951+952+953, (118) 



b47r CriT ^^^^-^ ^ ,^^-^ 
+N^AxN^AYNf'BYN!'Bx)d2{m^o/, m^o/, m.^^^ me,') 

--m^c. m^o^A^2Ax^iAV^iV^iBxC^o("^x'A^ m^o^^ mex^ mev')}, (119) 

^52 = -T7^^£[E E^2"4xiV/'Bx{4(%)Ai.C2(m,/,m^<,^m^.,2) 

-■-"^^ A,B=1X=1 



-2m^o m^o (K-o )^BCo(me^^ m^o 2, m^o 2)} 
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+ T. H N^AxNl'lYiXeJxYC2{m^o/,m,J,m,^% (120) 

A=l X,Y=1 

9h = 933, (121) 

/2 2 3 

^^'^ A,B=1X=1 

-2^x:,^x,(^x«)-4BCo(m^^',m;^-',m^^')}, (123) 

^53 = ^33- (124) 

Various mixing matrices and Z coupling constants which appear in the above for- 
mulas are given as follows: 



(%)aB = --{iONL)A3{ONLTB3- iONL)A^{ONLTB4}, (125) 



0^xI)aB = 7;{iONR)A3{ONRTB3 ~ {O N r) A4{0 N r)*B4} , (126) 



iy<,-)AB = --{OcL)A2{OcLyB2, (127) 

{y^^)AB = -liOcR)A2{OcRrB2, (128) 



{X.JxY = -j:iUe)xkm*Y,, (129) 

k=l 

iX,^)xY = E(f^e)xfc+3(f/e)yfc+3- (130) 



k=l 
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Zl = (-i + sin^M, (131) 

Z% = sin^ 9w (132) 
The photon-penguin couphng constant is written as follows: 

An = Al + A^j„ (133) 



^R ~ nKci^2r^_ 51 Yl —2-{a^2AxNiAxbi 



2567r2Gi7' j^^ ^^ ml^ 6 m- 



2 



+N,\xN,^Ix^^bU^)}: (134) 



a/2p 2^11 "^-- 



A=1X=1 ""yx " ""i'x 

m— m^— 
+C^AxC?lx^^b^2{^)}- (135) 



The other coupling constants can be obtained by simply exchanging the suffix of 
above formulas: 



52 = 9i{L^R), (136) 

94 = 93{L^R), (137) 

g, = 95{L^R), (138) 

Al = Ar{L^R). (139) 

B.3 Mass Functions 

The mass functions used in the effective coupling constants of the ^^ -^ e+7 and 
/i"*" -^ e^e^e" processes are defined as follows: 
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K{x) 



bUx) 

hl{x) 
bUx) 



hl{x) 



1 



2(1 -x)4 

{i-xY^ 

{l-xf^ 

1 

2(l-a;)4 

1 

2(l-a;)4 

1 

2ri-a;)3 



(2 - 9a; + ISx^ - llx^ + 6a;^ ln(a;)), 

1 - 6x + 3x^ + 2x^ - Gx^ ln(x)), 

1 -x^ + 2a;ln(x)), 

(-16 + 45a; - "iQx^ + Ix^ + 6(3a; - 2) ln(a;)), 

(2 + 3a; — 6a;^ + x^ + 6a; ln(x)), 

(-3 + 4a;-a;2-21n(a;)), 



(140) 
(141) 
(142) 
(143) 
(144) 
(145) 



co{x,y,z) 
C2{x,y,z) 



xln(x) 



y\n{y) 



zln{z) 



{y-x){z-x) {x-y){z-y) {x - z){y - z)' 



1.3 x^\n{x) 



y'^My) 



z^ 111(2;) 



(146) 



4 2 {y-x){z-x) {x-y){z-y) {x - z){y - z) ' 



(147) 



doix,y,z,w) 



xln(x) 



+ 



y\n{y) 



{y - x){z - x){w - x) {x-y){z~y){w-y) 
z\n{z) w\n{w) 



+ 



{x — z){y — z){w — z) (x — w)(|/ — w){z — w)' 



:i48) 



d2{x,y,z,w) 



x"^ ln(x) 



t/2 ln(t/) 



4 (y — x) (2; — x) (w — x) (x — y) {z — y){w — y) 



z'^\n{z) 



+ 



w^ln(w) 



(x — z)(|/ — z){w — z) (x — w)(|/ — w)(z — w) 



}. (149) 
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C Neutron EDM 

We discuss QCD correction in the calculation of the neutron EDM [^,|2ll]. The 



neutron EDM are calculated by the following effective Lagrangian. 

^cfT = E Cf (/i)Of (/i) + E Cf (/x)Of (/i) + C«(^)0^(/i), (150) 



q \r~;- q \r~ / • z^ ^q 



where O^, Og, O^ correspond to the quark electric dipole, chromomagnetic dipole, 
and gluonic Weinberg's operators, respectively, which are given by 

Of = -'-qa.^^.qF^'^, (151) 

O^, = -'^qcr.ul.T^qG^^'', (152) 

O^ = -h"'"(^''"^'GlpGl^Gl''. (153) 

Here, e°^^^ = 1, and /"'"^ is the structure constant of the SU(3) group. 

In SUSY models, we can obtain the Wilson coefficients at the electroweak 
scale by evaluating 1-loop diagrams. Cf is induced by the photon-penguin diagram. 
C^^^ -' and G^^^ -* are induced by the photon- and gluino-penguin diagrams. There 
are three types of SUSY contribution, chargino-squark loop, neutralino-squark loop 
and gluino-squark loop diagrams. The gluonic Weinberg's operator is induced at 
a 2-loop level and the diagram involving the stop and the gluino gives dominant 
contribution. These contributions are listed in Ref.pl . 



We can take into account a QCD correction from the electroweak scale to a 
hadronic scale (1 GeV), by using the following renormalization group equations for 
the Wilson coefficients. 

i^ = ^-,-C(,), (154) 

an Air 

where C = (Cf, Gg, G^)'^ and the anomalous dimension matrix 7^^ is written by 

/ 8/3 

7 = 32eg/(3^,) (-29 + 2Nf)/^ | . (155) 

V Qrriq 2A^/ + 3 



33 



Here, Nf is a number of the quark flavor and Q denotes the electro-magnetic charge 
of the quark in unit e (e > 0). The RGEs can be solved analytically as follows: 



C^f(/i) 



T]'^ 



C^M + 8eQ[l-v 



C^if^) 



72eQmq{fio) 

7 + 2Nf 

-29+2]Vx 



1 — r/' 



V + 



2iVf + 5 






fi's(/^o) 



(156) 



V 



c::{fi, 



9 



14+4NJ-- ' 



7 + 2N. 



l_^33-2iv, Wn,(/io)C^(/xo) 



/ 



C^(/x) 



9+6Nf 
33-2iVf /^G 



C^^(/io), 



(157) 
(158) 



where rj = 5f,(/io)/5's(/w). 

We solve RGE from rriw to ttij,, ruf, to mc and r/ic to the 1 GeV scale. When 
the heavy quarks (c, b) decouple at their mass threshold, C^ is induced through the 
chromo-electric dipole moment of the heavy quarks. Difference C"^ below and above 
the threshold is given by [^ 



C^(m,) 



below 



-C^im 



q) above 



"•'"'> <(»,). 



■^7rmq{mq 



(159) 



Taking into account the QCD and threshold corrections, we obtain the effective 
Lagrangian at the hadronic scale. It is then straightforward to evaluate the effective 
£ at 1 GeV scale from rriw scale. 

The neutron EDM (dn) is given by the Wilson coefficients at a hadronic scale 
as follows: 

(160) 
(161) 

(162) 



dn = 


= d^ + d^ + d^, 


d^n - 


- li^c^-c-), 


d'n - 


^ ^£K-^"^)' 


d^ - 


eM G 
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where M is a chiral symmetry breaking parameter, which is estimated as 1.19 GeV. 
In the above we use non-relativistic quark model for d^ and naive dimensional 
analysis for d'^ and d'^. 
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Figure Captions: 
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FIG. 1: Kinematics of the /i^ -^ e^e^e" decay in the center-of-mass system of 
muon. The plane I represents the decay plane on which the momentum vectors pi, 
P2, ps lie, where pi and p2 are momenta of two e+'s and p^ is momentum of e^ 
respectively. The plane II is the plane which the muon polarization vector P and p^ 
make. 
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0.5 



X 



FIG. 2: Kineinatical Region of the /i^ -^ e^e^e~ decay in the center-of-mass system 
of muon. Xi (X2) represents a larger (smaller) energy of decay positrons normalized 
by ^. A light shaded region is allowed. 
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FIG. 3: The observables in the SU(5) model without the SUSY CP violating phases 
in the me^-|Ao| plane. We fix the SUSY parameters as tan/? = 3, M2 = 150 GeV 
and /i > and the top quark mass as 175 GeV. (a) Branching ratio for fi'^ -^ 
e~^'~f normalized by |A^p= |(Vr)23(Vr)^3P. (b) Branching ratio for fi^ -^ e+e+e^ 
normalized by \\t\'^. (c) The ratio of two branching fractions ^(^J'l ■ (d) The P- 
odd asymmetry for yU+ -^ 6^7. (e) The P-odd asymmetries Ap^ for /i+ — > e^e^e". 
(f) The P-odd asymmetries Ap^ for /x"*" — > e~^e^e~ . The cut-off parameter 6 is taken 
to be 0.02. 
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FIG. 4: ^^0 and d^ dependences on the EDMs and A^. We take a specific set of 
SUSY parameters tan/3 = 3, M2 = 300 GeV, rrag^ = 650GeV and |y4o| = 1 in tlie 
parameter region — vr < Qa^ < ^r and — O.OSvr < 9^ < O.OSvr. Dark sfiaded regions 
are excluded by tlie EDM experiments 
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FIG. 5: The obsevables in the SU(5) model with the SUSY CP violating phases in 
the me^-|2lo| plane. We fix the SUSY parameters as tan/5 = 3, M2 = 300 GeV, 
6ao = f and 6^ = Q and the top quark mass as 175 GeV. (a)-(f) are same as Fig. |^. 
(g) The T-odd asymmetry for /i+ — *■ e+e+e". The cut-off parameter 5 is taken to be 
0.02. The experimental bounds from the electron, neutron and Hg EDMs are also 
shown in each figure. The left upper solid line corresponds to the electron EDM, 
the right upper solid line to the neutron EDM and the right lower solid line to the 
Hg EDM. The lower side of each bound is allowed by these experiments. The upper 
side of the bold line is excluded by the EDM bounds even if we allow 9^i_ taking 
slightly different value from 0. 
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FIG. 6: The obserbables in the SU(5) model for tan/3 = 10 in Tng^-IAol plane. Other 
parameters are same as in Fig. |^. 
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FIG. 7: The observables in the SO(IO) model without the SUSY CP violating phase 
in me^-|ylo| plane. The input parameters are same as in Fig. ^. The upper right 
corner bounded by the bold line is excluded by phenomenological constraints and 
the small rrigj^ region bounded by the left bold line is not allowed in the minimal 
SUGRA model. 
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FIG. 8: The obsevables in the SO(IO) model with the SUSY CP violating phase in 
r?2e^-|Ao| plane. The input parameters are same as in Fig. ^ The small mg^ region 
bounded by the left bold line is not allowed in the minimal SUGRA model. The 
upper right bold line shows the same bound as in Fig. H. 
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FIG. 9: The differential branching ratio and asymmetries for the /i"*" -^ e^e^e^ 
process in the SU(5) model as a function of x\ which is a larger energy of decay 
positrons (^^). We fix the SUSY parameters as tan/5 = 3, M^ = 300 GeV, mg^ = 
TOOGeV, |y4o| = 0.5, Oa^ = f and 6'^ = 0. (a) The differential branching ratio for 
the /i+ -^ e'^e+e^ normalized by the total branching ratio, (b) The differential P- 
odd asymmetry ap-^ . (c) The differential P-odd asymmetry ap^ . (d) The differential 
T-odd asymmetry a^. 
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FIG. 10: The differential branching ratio and asymmetries for the /i"^ —>■ e+e^e^ 
process in the SO(IO) model as a function of Xi. The input parameters are same as 
in Fig. ^ 
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